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The Z^-invariant chiral Potts model is considered as a perturbation of a Z^ 
conformal field theory In the self-dual case the renormalization group equations 
_. . become simple, and yield critical exponents and anisotropic scaling which agree 

^ . with exact results for the super-integrable lattice models. The continuum theory 



is shown to possess an infinite number of conserved charges on the self-dual line, 



c^ ■ which remain conserved when the theory is perturbed by the energy operator. 
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The Zf^ chiral Potts models, first introduced as models of systems exhibiting 
commensurate-incommensurate and melting transitions [1,2], have more recently 
become of interest as exactly solvable lattice models whose Boltzmann weights 
satisfy the star-triangle relation, yet do not enjoy the so-called 'difference property' 
possessed by previously studied such models [3]. The chiral Potts models are also 
unusual in that they are intrinsically anisotropic in their critical properties [4,5], 
unlike most other such lattice critical systems in which rotational symmetry is 
recovered in the scaling limit, so that they are described by a conformal field 
theory. 

A special case of the chiral Potts model, the so-called superintegrable case, 
has proved to be especially tractable [2,5,6,7]. It becomes critical at the self-dual 
point, and the critical exponents describing the specific heat singularity and the 
anisotropic scaling of the correlation lengths have been obtained exactly: a. = 
1 — 2/A^, v\ = 1, and U2 = 2/N [4,5]. In addition, the presumably exact conjec- 
tures (3j = j{N — j)/2N have been made [2,4,8] for the exponents governing the 
order parameters (e*-^ ) ~ {—t)^^. Interestingly enough, the values found for the 
'thermodynamic' exponents a and Pj agree with those of the isotropic Z^ models 
[9] obtained when the chirality parameter is set to zero, whose continuum limit at 
criticality is believed to correspond to the Z^y conformal field theories first studied 
by Fateev and Zamolodchikov [10]. From the renormalization group point of view, 
this result is puzzling, since, as will be discussed below, the chiral perturbation of 
the isotropic model appears to be relevant, and therefore should be expected to 
modify the critical exponents in an essential way. 

In this paper, we give a simple explanation of this apparent paradox, showing 
why the thermodynamic exponents are not changed. In addition, we give a simple 
derivation of the anisotropic scaling exponents ui and z/2- Our methods are strictly 

2 



valid only for weak chirality, close to the isotropic point. However, they suggest that 
the exponents are in fact independent of the magnitude of the chirality parameter, 
to all orders in perturbation theory. The main non-universal feature, which appears 
to vary continuously along the critical line, is the orientation of the principal axes 
with respect to which anisotropic scaling occurs. In addition, we show that there 
is a two-dimensional parameter space in which these models are integrable in the 
continuum sense, that is, they possess an infinite number of conserved charges. 

In the simplest case, A^ = 3, the reduced lattice Hamiltonian for this model on 
a square lattice with sites labelled by pairs of integers (a;,y) is 

5 = - ^ {Kx cos{AxOx,y - 5x) + Ky cos{Ay0x,y - 5y)) , (1) 

x,y 

where the 9x,y are angles which are integer multiples of 27r/3, and A^ (n = x, y) is 
the lattice difference operator. The original version of this model, due to Ostlund 
[1], corresponds to dy = 0. The self-dual case, of interest here, occurs for dy = ±idx- 
The superintegrable point then corresponds to 5a; = ^. When dx = Sy = 0, this 
model is the ordinary 3-state Potts model, whose continuum limit at the critical 
point is described by the Z3-invariant conformal field theory with central charge c = 
|. When the chiral perturbation is turned on, to which operators in the continuum 
theory does it couple? On the lattice, S^ couples to sin A^9. This transforms under 
90° lattice rotations like a vector, and so should, in the continuum limit where full 
rotational invariance is recovered, transform as a combination of operators with 
spins ±(4?i + 1), where n = 0, 1, . . .. The most relevant such operators of the Z^ 
conformal field theory, which are not total derivatives, have spin ±1, and (labelled 
by their scaling dimensions {h,h)) are $/? 2\ and $/2 i\. Note that they have 

V 5 ' 5 '' V 5 ' 5 '' 

overall scaling dimension x = | < 2. This means that they are relevant in the 
renormalization group sense. Considerations of rotational symmetry then imply 
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that S = Sx + iSy couples to $/7 2-> and d = dx ~ iSy to $(-2 7\. 

We now consider the duahty properties of these operators. The duahty transfor- 
mation involves three steps: first the Boltzmann weights on each link are Fourier 
transformed with respect to A^9, introducing new integer-valued conjugate link 
variables n^, then the sums over the original variables 9x,y are carried out, which 
give constraints of the form A^n^ = (mod 3) at each lattice site. These are then 
solved by writing n^ oc e^i,9iy, where the 9iy are angles defined on the dual lattice, 
and e^// is the totally antisymmetric symbol. Since, in the Fourier transform, n^ 
couples linearly to z5^, it follows that 9i, will couple to ieu^d^. Keeping track of 
the factors more carefully, we find that under duality, 

V -* i^iiv^v , (2) 

so that 6 ^> 6 and 6 -^ —d. The model with 5 = 0, 5 7^ is thus self- 
dual. Equivalently we may say that, under duality, $^7 2^ -^ $/7 2\, while* 
^(2. 1) ^ ^^(2. 1\- This is consistent with the fusion rules of these operators, 
$/7 2\ ■ ^(2 7^ ~ $/2 2\ -(- ■ ■ ■, since the leading thermal operator $^2 2\ is odd 
under duality. 

Thus we are led to consider the field theory described by the action 

S = So + 6 / $.7 2^(fix + 6 $.2 Z^(Px + T $/2 2-.(Px , (3) 

J VS'S-* J VS'S-* J ^S'S/ 

where 5*0 is the action of the Z^ conformal field theory. The self-dual subspace, 
with which we shall be mostly concerned, then corresponds to 5 = r = 0. We are 
therefore considering a conformal field theory perturbed by an irreducible spin one 
operator. This will turn out to have greatly simplifying consequences. 



* 



This apparent lack of symmetry between $/7 2n, and $/2 7^ can be traced to a particular 

choice for e^jy. In fact, there are two possible duality transformations corresponding to 
each choice of sign. Under the other kind, it is $/2 7\ which is even. 



Let us consider the renormalization group equation for the dimensionless cou- 
phng g = 5a /, reflecting the flow of this couphng under changes in the microscopic 
lattice cut-off a. It is conventional to work within a scheme where the lattice cut- 
off is replaced by one which respects the rotational symmetry of the theory. In 
particular, we may consider an expansion of the partition function in terms of 
the connected correlation functions of $/7 2^, with an ultraviolet regulator which 

V 5 ' 5 '' 

prohibits any two of the integration arguments approaching more closely than a. 
Within such a scheme, the coefficients of the renormalization group /3-function to 
any finite order in perturbation theory are given by integrals over correlation func- 
tions evaluated in the unperturbed theory, which has rotational symmetry. This 
has the consequence that, to all orders in perturbation theory, the renormalization 
group equation for g is given by its first term only: 

9 = l9 , (4) 

because any higher order terms g^ with n > 1 would not transform correctly under 
rotations. Explicit calculations at low orders confirm this result. In the same way, 
if we consider any scalar perturbation of the action S, for example the thermal 
operator, then its renormalization group equation can contain no dependence on 
(7, so that 

f = |r + 0(r3) . (5) 

The same is true for all other such scalar quantities which could enter the free 
energy. 

The renormalization group equations, within such a scheme, are therefore triv- 
ial. There are only two possible fixed points, g = (the ordinary 3-state Potts 
model) and (7 ^ oo, to which all theories with g ^ fiow. Note that this fixed 
point lies at infinity in a rotationally invariant regularization scheme, but in others, 
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for example a lattice scheme where terms 0{g ~^ ) may appear on the right hand 
side of Eq. (4), it may move into a finite value. The second equation Eq. (5) implies 
that the renormalization group eigenvalues of all scalar operators are the same at 
g = and g ^ oo. This has the consequence that the exponents governing the 
singular behaviour of the free energy as a function of temperature, magnetic field, 
and so on are the same. It is consistent with the observation that they agree also 
at the superintegrable point. 

Next, we discuss the anisotropic scaling of correlation functions. This is a 
consequence of the appearance of a novel kind of space-time symmetry, which is 
most easily seen through the properties of the stress tensor. The argument is a 
simple generalization of that given by Zamolodchikov [11] for perturbations of a 
conformal field theory by a scalar operator. At the conformal point 6 = there are 
two non-zero components T = Tzz and T = T^z, where dzT = dzT = 0. The other 
components Tzz and Tzz vanish. When 5 7^ 0, this is no longer true. A correlation 
function {T{z) . . .) is modified to first order in 5 by a term 

-5j{T{z)^^7_2^{zi,zi) )d\ . (6) 

Using the operator product expansion 

7 i_ I 

T{z)-^(7_2,{zi,zi) = - — ^— 2*(Z2^(^,^) + ^-dz^(7 2,{z,z) + --- , (7) 

and the fact that dz{z — zi)^ = ^6^ '{z — zi), it follows that 

c>,-T =1(1- I) 5 $(7^2) +0(52) . (8) 

Now we should expect to be able to express all the terms on the right hand side 
in terms of the renormalized coupling constant and renormalized operators of the 
conformal field theory. But Eq. (4) implies that 6 needs no renormalization (as 
may be checked, there are no new ultraviolet divergences in correlation functions), 
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and a simple dimensional argument shows that the operator multiplying the term 
of order 5" would have scaling dimensions (2 + 2n/5, 1 — 3n/5). Since these must 
all be non-negative, only the n = 1 term is in fact present. It then follows from 
the conservation equation dzT + dzTzz = that 



5(1-1)^^(12) : (9) 

^ V 5 7 5 7 



and similarly that 



2 



-f(l-|)5$(y) . (10) 



The fact that these are not equal implies that rotational symmetry is broken, as 
expected. Instead we see that 

Tzz = -\Tz-z . (11) 

This implies the existence of an unusual symmetry, under transformations which 
are mixtures of rotations and dilatations 

z^e^'^z , z^e^'^z , (12) 

where u is arbitrary. However, we see immediately that these transformations do 
not preserve the reality of z + z and i{z — z), and therefore do not appear to be 
physical. This requires careful discussion. Since 5y is pure imaginary, the Boltz- 
mann weights are not necessarily real and positive. This means that there is no 
straightforward interpretation as a statistical model. Several of the usual prop- 
erties of such systems, such as insensitivity of bulk thermodynamic quantities to 
boundary conditions, can fail in this case, and it appears that the usual arguments 
supporting the existence of a continuum limit may not go through either. How- 
ever, as pointed out by several authors, although the transfer matrix corresponding 
to Eq. (1) is not real, it is hermitian, as is the quantum hamiltonian which arises 
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in the limiting case. Therefore such a transfer matrix, or quantum hamiltonian, 
should define a sensible quantum theory in 1 + 1 dimensions with the possibility 
of a continuum limit as a quantum field theory.* Within this picture the symme- 
try Eq. (12) is now physical, since in real time z and z get replaced by light-cone 

co-ordinates x = t ± x, while Eq. (11) becomes T |- = ^T-| The generator of 

this symmetry is ^D + \^-> where D is the generator of dilatations, and M that 
of Lorentz boosts. 

This has important consequences for the spectrum. If the hamiltonian and 
momentum operators are denoted by H and P, and we introduce their light-cone 
combinations P = H ± P^ then it follows from the commutation relations of 
if, P, D and M that if |p^,p^) is an eigenstate, so is |e^ p+,e^ p^). The 
density of states is constant along curves p^ = nd^ ' {p'^) ' \ where k is constant. 
Within perturbation theory in 5, only states with positive p'^ and p^, and therefore 
positive K, may arise, since the unperturbed correlation functions have support 
only inside the light-cone. However, this does not rule out the possibility of non- 
perturbative states with k < 0. For the spectrum of H to be bounded from below, 
there must be a minimum value of k, say ^niin- If /^min < 0, the ground state 
will then lie on this curve, and will have non-zero eigenvalue of P, implying a 
spontaneous breaking of translation symmetry. In fact, exact calculations at the 
superintegrable point [12], and numerical calculations for more general values of 
the chiral parameter [13], indicate that this in fact happens. As the critical line is 
approached, a whole sequence of first-order transitions takes place in which states 
with non-zero momentum cross the lowest energy state in the zero-momentum 



* These observations do not invalidate our previous assertions about the structure of the 
renormalization groupequations, which were made within the context of perturbation the- 
ory, which is well-defined in Euclidean space. However, it should also be possible to work 
entirely in real time, and require a perturbative regularization procedure which respects 
Lorentz invariance. 



sector. The results for the critical exponents, however, refer to quantities evaluated 
in the perturbative sector, whose lowest energy state we denote by |\I/o)- 

A two-point correlation function of an arbitrary scalar operator (j) can be written 
as an integral over a spectral density 

(vl/o|T[(/>(x+,x-)(/>(0,0)]|vl/o) = J pip+,p-)e'P-'^dp+dp- , (13) 

where 

Pip+.p-) = J2\('^0\<P\n)\^S{p+-p+)6{p;,-p-) (14) 

n 

which will have the scaling form 

p(p+,p-) = (pV)'"^-Vf^^^j , (15) 

where we have used the facts that (p transforms under dilatations according to its 
scaling dimension x^, and is a Lorentz scalar, and that |\l/o) is annihilated by P^ 
and P~ . 

All this is supposed to be valid on the self-dual line. When the thermal pertur- 
bation is added, the scaling function $ will also depend on the additional scaling 
variable r(p+p^)^ ' (where the exponent | reflects the value of the thermal eigen- 
value.) We may write this two- variable scaling form alternatively as 

p{p+,p-) = (pV)"^"V~(p+r,P"e+) , (16) 

where 

(17) 

We therefore identify the anisotropic correlation length exponents z/i = 1 and 

1^2 = 2/3. 

It will be noticed immediately that the principal axes with respect to which 
this anisotropic scaling takes place are x , not t and x as is found in studies of 
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the superintegrable case [2,4]. We now offer an explanation of this. The lattice 
theory differs from the continuum action Eq. (3) by terms which break rotational, 
or Lorentz, symmetry. These terms are irrelevant in the isotropic theory at the 
critical point, since they correspond to operators of spin 4 or higher. When such an 
isotropic critical theory is perturbed by a scalar operator, it is generally true that 
such operators remain irrelevant, at least in the scaling region. This is consistent 
with the observation that no new terms in their renormalization arise from the 
scalar perturbation. When a theory is perturbed by a non-scalar operator, however, 
this is no longer valid. Indeed, in the present case the coupling constant Qs of an 
operator of spin s can be renormalized at order g^ by the interaction, so that its 
renormalization group equation has the form 

gs = ys9s + const. 5f^ H . (18) 

A particular case of interest is the (zz) component of the stress tensor T whose 
coupling (72 will get a contribution at order g . In fact, this may be computed from 
the operator product expansion $/7 2\ ■ $^7 2\ ~ AT + ■ ■ ■ to give 

92 = ^9^ + --- . (19) 

The stress tensor is an example of a redundant operator. When added to the action, 
it does not modify the critical exponents, but instead its effect may be removed 
by an appropriate co-ordinate transformation. In our case, this is accomplished by 
sending* 

x+ ^ x+ + 7n'^g'^x~ . (20) 

There are, in addition, higher order effects in a which come from terms in Eq. (19) 
like g^4g J and so on. All of these will be non-universal, and cannot be calculated 



* The operator product expansion coefficient ^ comes from the scaling dimension ^ normal- 
ized by § = ■5- The extra factor of 27r in Eq. (20) arises from that in the conventional 
definition of the stress tensor. 
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within our continuum approach. Note that aU these additional operators do not 
affect the simple renormalization group equation Eq. (4) for g, since they are all 
have even spin. Although this shows that non-universal lattice effects may rotate 
the principal scaling axes away from the light-like directions x , it does not explain 
why the axes should be precisely parallel to t and x at the superintegrable point. 
Perhaps there is a hidden symmetry at this point which requires such an alignment. 

Although the arguments above have been presented for A^ = 3 for clarity, they 
generalize straightforwardly to arbitrary A^. In that case, in the isotropic Z^y 
model, the leading thermal operator has scaling dimensions (2/(A^ + 2), 2/(A^ + 2)), 
and the leading spin 1 operator ((A^ + 4)/(A^ + 2), 2/(A^ + 2)). We take this as 
coupling to the chiral parameter 5. All the previous arguments then go through, 

and in particular we find that T |- = jjT^ This leads to anisotropic correlation 

length exponents z/^ = 1, z/2 = 2/A^, in agreement with the exact results found at 
the superintegrable point [4,5]. 

Finally, we give a simple argument that these continuum models should be 
integrable. In the conformal field theory, there is an infinite number of conserved 
charges corresponding to holomorphic and antiholomorphic currents {T^^\z)^Q) 
and (0, T {z)) of arbitrarily high spin s. When the theory is perturbed by the 
thermal operator $2/iV+2 2/N+2^ ^^ ^^ ^^^^ known that an infinite subset of these 
remain conserved, with spins s following the exponents of of A^y-i modulo the 
Coxeter number [14]. This may be traced to the property that the residue at 
zi = z in the operator product expansion of T^^'{z) and ^2/N+2 2/N+2i^l^ ^l) ^^ ^ 
total derivative with respect to z [11]. Similarly, an infinite set of antiholomorphic 
currents remain conserved for an analogous reason. Now all of the Zj^ conformal 
field theories possess an extended VF-symmetry, generated by conserved currents 
W^^' (z) and W (z) [15] . With respect to this extended symmetry, the chiral fields 
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of dimension {{N + 4)/(A^ + 2), 2/{N + 2)) are not primary, but VF-descendents 
of the form W_-^^2/N+2 2/N+2 ^^^ some n. Since the holomorphic current W^^' 
commutes with all the antiholomorphic ones, it follows that all the antiholomor- 
phic currents which remain conserved under a purely thermal perturbation remain 
conserved when the chiral perturbation is also present. In general, however, the 
holomorphic currents are no longer conserved. It is unclear at this point what con- 
sequences may be drawn from this. The theory with both types of perturbation 
is expected to be massive, and therefore to possess a sensible 5'-matrix. However, 
in the absence of Lorentz invariance (which holds only when 5 = 0) , or the mixed 
symmetry discussed earlier (which holds only when r = 0), there is no simple way 
of determining how these conserved charges act on the asymptotic states. 

The manifolds of integrability are the planes 6 = and 6 = in the space 
parametrized by (r, 5, d). These appear to differ from the lattice case, where, close 
to the isotropic point, the exact solution manifold, expressed in our coordinates, 
has the form r ~ 66. It is possible that irrelevant couplings are responsible for this 
difference, or that continuum and lattice integrability are rather different proper- 
ties. 

In passing, we comment on the more physical case of of the uniaxial model with 
6x ^ 0, 6y = 0. In this case symmetry arguments are are far less powerful. The 
renormalization group equation for 6x is no longer trivial, and gets corrections at 
order 6^. These are very difficult to compute, involving as they do an integral over 
a 6-point function. However, they do suggest the possibility of a non-trivial fixed 
point for 6x 7^ 0. If this were to occur, it would correspond to a Lifshitz point, 
whose existence is still controversial in numerical studies of the A^ = 3 model [16]. 

In summary, we have given simple symmetry arguments for the values of the 
critical exponents of the Z]\[ chiral Potts models on the self-dual line. These argu- 
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merits are valid in the limit of weak chirality, but the results are independent of the 
its magnitude, and are consistent with those found at the superintegrable point. 
We have argued that the continuum theory makes sense only in Minkowski space, 
and found that it exhibits an unusual space-time symmetry, which allows for the 
observed spontaneous breaking of translational symmetry in the true ground state. 
Our analysis indicates that, while the anisotropic correlation length exponents are 
universal, the principal axes to which they refer are not, and that these should 
rotate away from the light-like directions as the chirality parameter 6 is increased. 
This may explain the inconclusive results obtained for mass-gap scaling in numer- 
ical work for small values of 6 [2] . Crossover theory suggests that the asymptotic 
regime should be reached only for reduced temperatures of the order of 5 ''^ ^^, 
which will be a very narrow region, especially for A^ = 3. However, even then, 
the exponent z/^ = 1 should be observed only if the gap to states with the correct 
non-zero value of the momentum is considered. In all other cases, the other ex- 
ponent U2 = 2/N should be observed. It would be interesting to see whether the 
predicted rotation of the scaling axes occurs in the vicinity of the superintegrable 
point, where numerical work should be more conclusive. 

The author thanks R. Baxter and B. McCoy for useful discussions on the lattice 
results, and the Isaac Newton Institute for Mathematical Sciences for its hospital- 
ity. This work was supported by the Isaac Newton Institute, the UK Science and 
Engineering Research Council, and by the US National Science Foundation Grant 
PHY 91-16964. 
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